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Abstract 

In the scenario with four quark generations, we perform a fit using flavor-physics data 



and determine the allowed values - preferred central values and errors - of all of the 
elements of the 4x4 quark mixing matrix. In addition to the direct measurements of 
some of the elements, we include in the fit the present measurements of several flavor 



changing observables in the K and B systems that have small hadronic uncertainties, and 
also consider the constraints from the vertex corrections to Z — > bb. The values taken for 



> 

en . 

the masses of the fourth-generation quarks are consistent with the measurements of the 

(N ■ 

oblique parameters and perturbativity of the Yukawa couplings. We find that \Vtb\ > 0.98 
at 3<t, so that a fourth generation cannot account for any large deviation of \Vtb\ from 
unity. The fit also indicates that all the new-physics parameters are consistent with zero, 
and the mixing of the fourth generation with the other three is constrained to be small: 
we obtain \V u y\ < 0.06, \V c y\ < 0.027, and \Vq/\ < 0.31 at 3d. Still, this does allow for 
the possibility of new-physics signals in Bd, B s and rare K decays. 

PACS numbers: 12.15.Hh, 14.65.Jk 
Keywords: CKM matrix, Fourth generation 



* Electronic address: 
^Electronic address: 
^Electronic address: 



alok@lps . umontreal . ca 
amol@ theory, tifr . res . in 



london@lps . umontreal. ca 



1 



I. INTRODUCTION 



There is no unequivocal theoretical argument which restricts the number of quark 
generations to three as in the standard model (SM). An additional fourth generation 
(SM4) is one of the simplest extensions of the SM, and retains all of its essential 
features: it obeys all the SM symmetries and does not introduce any new ones. At 
the same time, it can give rise to many new effects, some of which may be observable 
even at the current experiments |l| . Even though the fourth-generation quarks may 
be too heavy to have been produced at the pre-LHC colliders, they may still affect 
low-energy measurements through their mixing with the lighter quarks. The up-type 
quark if would contribute to b — > s and b — > d transitions at the 1-loop level, while 
the down-type quark b' would contribute similarly to c — > u and t — > c. 

The addition of a fourth generation to the SM leads to a 4 x 4 quark mixing 
matrix CKM4, which is an extension of the Cabibbo-Kobayashi-Maskawa (CKM) 
quark mixing matrix in the SM. The parametrization of this unitary matrix requires 
six real parameters and three phases. The additional phases can lead to increased 
CP violation, and can provide a natural explanation for the deviations from the 
SMpredictions seen in some measurements of CP violation in the 5-meson system 
{2 7]. A heavy fourth generation can play a crucial role in the dynamical generation 
of the electroweak (EW) symmetry breaking [8[]. Also, the large Yukawa couplings 
of the fourth generation quarks, together with the possible large phases, can help 
efficient EW baryogenesis 9]. 

The EW precision measurements of the oblique parameters S and T imply strong 
correlations between the masses of the fourth-generation quarks 10|, yjj . The pa- 



rameter space of fourth-generation masses with minimal contributions to S and T, 



and in agreement with all experimental constraints, is 
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12| 



m t > > 400 GeV 



m t ,-m b , ~ ( 1 + ~ ™* ) x 50 GeV , (1) 



5 115 GeV 

where mt'imy and m# are the masses of t', b', and the Higgs boson H, respec- 
tively. On the other hand, the perturbativity of the Yukawa coupling implies 
that m t / ^ y/2z(v) ~ 600 GeV, where (v) is the vacuum expectation value of 



2 



the Higgs. Arguments based on the unitarity of partial S-wave scattering ampli- 
tudes for color-singlet, elastic, same-helicity t'-i' scattering at tree level restrict 



mf ^ ^/Ait/3(v) ~ 500 GeV 13|, |l4{. Thus, the fourth-generation quark masses 
are constrained to a narrow band, which increases the predictivity of the SM4. 

The quark-mass bounds above may be somewhat relaxed with the introduction 
of heavy fourth-generation leptons, which help in partially cancelling out the effect 
of the fourth generation on the S and T parameters. Even in the absence of any 
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quark-lepton cancellation, the EW precision measurements restrict 

m v - /TV ^ (30-60) GeV , 

where my and m v > are the masses of the fourth-generation charged lepton /' and 
neutrino z/, respectively. Thus even in the absence of any fine-tuned cancellations, 
there is a significant allowed range for the masses of the fourth-generation fermions, 
which is, in fact, not beyond the reach of the LHC The invisible decay width of 
the Z boson constrains the mass of the fourth-generation neutrino to be greater 
than 45 GeV. Though one would need a special mechanism leading to a massive 
fourth-generation neutrino and three ultralight SM neutrinos, phenomenologically 
this is perfectly allowed. 

In order to make concrete SM4 predictions, the first step is to determine the 
elements of CKM4. This involves not only fixing the values of the new parameters, 
but also re-evaluating those of the SM. This is because not all elements of the 
CKM matrix are measured directly. For example, the bounds on \V t d\ and \V ts \ are 
obtained from decays involving loops, and these are rather weak. And though |V#>| 
is measured in the tree-level decay t — > bW, its value is not that precise: the direct 
measurement at the Tevatron from single top production gives |T4&| = 0.88 ± 0.07 
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17fl . Now, \V t b\ = 1 is predicted in the SM to an accuracy of 10~ 3 . Although the 
Tevatron value is consistent with the SM prediction, it can also be as small as 0.67 
at 3cr. Thus, the values of the elements V tq (q = d, s, b) are not obtained through 
measurements. Rather, they are mainly determined using the unitarity of the 3x3 



CKM matrix 18]. However, the assumption of the unitarity of the 3x3 matrix is 
clearly invalid in the four-generation scenario, and relaxing it allows a much larger 
range of values for the elements \V tq \. For example, a large deviation of \V t b\ from 



unity is claimed to be possible in the SM4 [19M23 ] . 

We parametrize the CKM4 with 9 parameters, and perform a combined fit to 
these parameters using flavor-physics data. In addition to the direct measurements 
of the CKM4 matrix elements, the fit includes observables that have small hadronic 
uncertainties: (i) and from Z — > bb, (ii) ex from — > tttt, (iii) the branching 
ratio of K + — > 7T + uu, (iv) the mass differences in the Bd and B s systems, (v) 
the time-dependent CP asymmetry in Bd — > J/i/jKs, (vi) the measurement of the 
angle 7 of the unitarity triangle from tree-level decays, (viii) the branching ratios 
of B — y Xg'-f and B — > X c eu, and (ix) the branching ratio of B — > X s fi + fi~ in the 
high-g 2 and low-g 2 regions. We do not include the oblique parameters in the fit, but 
simply take the values of the fourth-generation quark masses to be consistent with 
the EW precision data. 
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'here have been several analyses of CKM4 in the past (e.g. see Refs. 5KZI, 
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24j]). However, they all have a number of deficiencies compared to the present 



work. They do not perform a fit. Instead, at best, they present scatter plots showing 
the allowed ranges of the CKM4 matrix elements (or correlations between various 
observables). Of course, these plots cannot quantify what the errors on the elements 
are, nor the confidence level of the ranges. This information can be obtained only 
by performing a true fit. Also, some of them do not include all clean observables 
which can be affected by the fourth generation. 

The paper is organized as follows. In Sec. HIl we define the Dighe-Kim 
parametrization of the CKM4 matrix. In Sec. IIII1 we present the observables which 
constrain the elements of CKM4, along with their experimental values. The results 
of the fit are presented in Sec. [TV] We conclude in Sec. [V] with a discussion of the 
results. 



II. CKM4 MATRIX: DIGHE-KIM PARAMETRIZATION 



The CKM matrix in the SM is a 3 x 3 unitary matrix: 



CKM3 



^ V u d V us V u b 
Vcd Vcs V c b 

v td v ts v tb 



(2) 



In the SM4, the CKM4 matrix is 4 x 4, and can be written as 



V ( 



CKM4 



(3) 



/ Vud Vus Vub V ub > \ 

V c d V cs V cb V cb i 
V t d V ts V t b V t y 
\ Vfd Vfs Vfb Vfb' J 

The above matrix can be described, with appropriate choices for the quark phases, 
in terms of 6 real quantities and 3 phases. 

In this paper, we use the Dighe-Kim (DK) parametrization of the CKM4 matrix 
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261 ] . This allows us to treat the effects of the fourth generation perturbatively and 



(4) 



explore the complete parameter space available. The DK parametrization defines 

V U8 = X, Vet = AX 2 , V ub = AX 3 Ce~ i5 - b 

V uV ee pX 3 e~ iS ^ , V cb > ee qX 2 e- i5 ^> , V w = rX , 

where A is the sine of the Cabibbo angle, so that the CKM4 matrix takes the form 

( # A AX 3 Ce~ i5ub P X 3 e- iS ^' \ 

# # AX 2 qX 2 e- iS ^' 

# # # rX 
V # # # # ) 

The elements denoted by can be determined uniquely from the unitarity condi- 
tion V r ( t KM4 VcKM4 = I- They can be calculated in the form of an expansion in powers 
of A such that each element is accurate up to a multiplicative factor of [1 + 0(X 3 )}. 



CKM4 



(5) 



The matrix elements V ud , V c d and V cs retain their SM values: 



V. 



ud = i - y + 0(X 4 ) , V cd = -A + 0(X 5 ) 



Vc, 



X' 



+ 0{X 4 



(6) 



whereas the values of the matrix elements Vtd, V ts and V tb are modified due to the 
presence of the additional quark generation: 

V td = AX 3 (1 - Ce iKb ) + rA 4 (qe i5 ^ - pe iS -»') 

+ ^A 5 (-r 2 + (C + Cr 2 )e iS " b ) + C(A 6 ) , 



Vts 

Vtb 



-AX 2 - qrX 3 e l ^ b ' + ^A 4 (l + r 2 - 2Ce^ b ) + C(A £ 



2 \2 



r z X 



+ 0(A 4 ) . 



(7) 



In the limit p = q = r = 0, only the elements present in the 3x3 CKM matrix 
retain nontrivial values, and the above expansion corresponds to the Wolfenstein 
parametrization j^] with C = \J p 2 + rf and 5 ub = tan" 1 (rj / p). 
The remaining new CKM4 matrix elements are: 

% d = A 3 (qe l5 c»' - pe iS ^') + ArX A (l + Ce l ^ b ) 

\5 

+ — (pe t8 ^' - qr 2 e iS ^' + pr 2 e iS ^') + £>(A 6 ) , 
V t , s = qX 2 e i5 ">' + ArX 3 

+ A 4 ^-pe i5 "»' + |e^' + ^e i5 ^ + £>(A 5 ) , 

% b = -rX + 0(X 4 ) , 

%v = 1 - ^ + 0(A 4 ) . (8) 



III. CONSTRAINTS ON THE CKM4 MATRIX ELEMENTS 



In order to obtain constraints on the CKM4 matrix elements, we perform ay 
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fit for all nine CKM4 parameters using the CERN minimization code MINUIT 28] . 
The fit is carried out for m t i = 400 GeV and 600 GeV. The b' mass is fixed by 
the relation m t i — my = 55 GeV [see Eq. (pQ)]. We include both experimental errors 
and theoretical uncertainties in the fit. In the following subsections, we discuss the 
various observables used as constraints, and give their experimental values. 

A. Direct Measurements of the CKM Elements 

The values of CKM elements obtained from the measurement of the tree-level 
weak decays are independent of the number of generations. Hence they apply to the 
3x3 and 4x4 matrices. The elements \V ur f\, IKJ, iKhl, \V rr ]\, \V r A and \V r h\ have 



all been directly measured. We use the following measurements 18j to constrain the 
CKM4 parameters: 

\V ud \ = (0.97418 ± 0.00027) , \V cd \ = (0.23 ± 0.011) , 

\V US \ = (0.2255 ± 0.0019) , \V CS \ = (1.04 ± 0.06) , 
\V ub \ = (3.93 ± 0.36) x 10~ 3 , \V cb \ = (41.2 ± 1.1) x 10~ 3 . (9) 
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B. Unitarity of the CKM4 Matrix 

Constraints on the CKM4 matrix elements can be obtained by using the unitarity 
of the CKM4 matrix. Through a variety of independent measurements, the SM 3x3 
submatrix has been found to be approximately unitary. We therefore expect all the 
CKM4 matrix elements which involve both the fourth-generation and light quarks 
to be relatively small. 

Using the measurements of \ V u d\, \ V US \ and \ V u b\, the first row of the CKM4 matrix 
gives 

iKbf = i - (IK/ + \v us \ 2 + \V ub \ 2 ) = 0.0001 ± 0.0011 . (10) 
Using the measurements of \V c d\, \V CS \ and \V c b\, the second row gives 

\V cb ,\ 2 = 1 - (\V cd \ 2 + \V CS \ 2 + \V cb \ 2 ) = -0.136 ± 0.125 . (11) 
Similarly, from the first column of CKM4, we have 

\V td \ 2 + \V t , d \ 2 = 1 - {\V ud \ 2 + \V cd \ 2 ) = -0.002 ± 0.005 . (12) 
Finally, the second column of CKM4 implies 

\V ts \ 2 + \V t , s \ 2 = 1 - {\V US \ 2 + \V CS \ 2 ) = -0.134 ± 0.125 . (13) 



Vertex Corrections to Z — > bb 



Including the QCD and QED corrections, the decay rate for Z — > bb is given by 



29] 



nZ -+ qq) = - . 2 a ™ Z 2R (\a q \ 2 + \v g \ 2 ) (l + if) 
16 sin V w cos^ 6w 

x (1 + Pqbd) (1 + Sqcd) (1 + ^) (1 + 5 q tQCD ) (1 + h) ■ (14) 

Here a q = 2J| and v q = ^2J| — A\Q q \ sin 2 Q w j are the axial and vector coupling 
constants, respectively. 

The S terms are corrections due to various higher-order loops: 

• 5 q ^ contains small electroweak corrections not absorbed in sin 2 6w- Their 
effect is at most at the 0.5% level. 
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Sq ED represents small final-state QED corrections that depend on the charge 
of final fermion. It is very small: 0.2% for the charged leptons, 0.8% for the 
w-type quarks and 0.02% for the d-type quarks 



Sqcd includes the QCD corrections common to all quarks; it is given by [29 1 

2 



>QCD 



a, 

— + 1.41 

7T V 7T 



(a. 



(15) 



where a s is the QCD coupling constant taken at the m z scale: a s = a s (m z ] 
0.12. 



5^ contains the kinematical effects of the external fermion masses, including 
some mass-dependent QCD radiative corrections. It is only important for the 
b quark (0.5%), and to a lesser extent for the r lepton (0.2%) and the c quark 

(0.05%) yiaoj. 



The correction 5^q CD consists of QCD contributions to the axia 
decay and originates from doublets with large mass splitting 
presence of the fourth generation, it is given by [32| 



29 



part of the 



3l|. In the 



$tQCD 



V q + a \ V 7F 



a. 



at /(/it) + at' f{lh>) + a>v f{Hb> 



where 



/(w)Klog( ^_,083 + !f5 + ^, 



(16) 



(17) 



with /i| = Am 2 Jm z . 

5b is non-zero only for q = b and is due to the Zbb vertex loop corrections. In 
the presence of the fourth generation, it is given by 32] 



8 h w 10 



-2 



2m| 



^r + 0.2 \V tb \ 2 + 



In order to isolate the large mass dependences appearing in the Zbb vertex 5b, 
one takes the following ratio 29|: 



R 



bb 



t 2 1 1 \-i 



(19) 



S 



where 



R s 
R c 



T(Z 66) 
T(Z -> ss) 
T(Z -> 66) 
r(Z -> cc) 

T(Z ->• 66) 

rrz ->• mm) 



0.9949 (l + , 

0.9960 1 1 + 1| ( 1 + (1 + 4), 
(1 + «2) (1 + V Q cd) { } 

i 0.9955-^ y~) *9££>h + s b ) 

(1 + «9 (1 + 8 c tQCD ) K b} 



(20) 



Using Eqs. (jlBjl-flSOjl. we get (for m t > = 400-600 GeV) 

-i 



bb 



i + 



3.584 
(1 + S b ) 



The data give 



33) 



R, 



bb 



0.216 ±0.001 



(21) 



(22) 



which, through Eq. ({TBI) , determines a linear combination of \Vtb\ 2 and |Vf/{,| 2 . This 
constrains the combination rX of the CKM4 parameters. 

We also consider constraints from the forward-backward (FB) asymmetry in Z — > 
bb. The Z — > bb interaction Lagrangian is 



C 



9 



h tl [9bLPL + 9bRP R )bZ i 



cos9 w 

where Pl(r) are the chirality projection operators, and 

1 1 



(23) 



9bL 
9bR 



§ + ^sin 2 6w + 6gl L + 6gt L , 
sin 2 6 W + 5g t bR + 5g t bR . 



(24) 
(25) 



Here the S's represent the radiative corrections due to the t and t' quarks. The FB 
asymmetry in Z — > bb allows us to determine the asymmetry parameter 1 



A h 



2 2 
9bL ~ 9bR 

9l L + 9m 



(26) 



1 The measured FB asymmetry A°p B is related to the asymmetry parameter via A°p B ss 
(3/ A)A e Ai,, where A e is the corresponding asymmetry parameter for the electron 34]. We only 
consider the parameter Ah since A e , and hence A°p B itself, would involve contribution from the 
fourth generation lepton sector, while we would like to restrict ourselves to the quark sector in 
this paper. 
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Within both the SM and SM4, only the cim, terms receive corrections proportional 



to m 2 t , at the loop level. We have 



22 



35- 



The data give [33] 
which, through Eq. 



bL 



bL 

9m = 0, 
9bR = 0. 





a 




m 2 


v tb \ 


16n sin 2 


6w cos 2 


9w 


m\ 




a 




m 2 , 


Vvi, 


167r sin 2 


9 W cos 2 


9 W 


m\ 



A h = 0.923 ± 0.020 



(27) 

(28) 
(29) 
(30) 

(31) 



constrains the combination rX of the CKM4 parameters. 



D. The K system 

Here we present observables in various K decays with the addition of a fourth 
generation. 



1. Indirect CP violation in Kl — > un 



Indirect CP violation in Ki — > tttt is described by the parameter ex, given by 

k e i( ^ £ 

£k = T^rhr, ; Im(M^) . (32) 



V2(AM K ) exp 

(AMjfjexp is the Kj-Kg mass difference. The parameters e = (43.51 ± 0.05)° and 



K e = 0.92 ± 0.02 39j include an additional effect from Im(v4 ), where A = A^K — )■ 
(7T7r)/ = o). M^ 2 is the off-diagonal element in the dispersive part of the amplitude for 
K°-K° mixing: 

= 2Z • (33) 



The calculation of Mf 2 in the SM4 gives [40 1 
'G 2 M 2 



M- 



K 
12 



f^f ) m K B K f K [r ]c (V; d V cs ) 2 S(x c ) + 2 Vct (V^^ 

+ Vt (v; d v ts ) 2 s(x t ) + 2 Vct iv^v cs )(v; d v t , s )s(x c ,x t ,) 

+ 2 Vtt iV; d V ts )(V; d V t , s )S(x t ,x t ,) + r ]t ,{V t 1 d V tls ) 2 S{x t ,)] . (34) 
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For the decay constant and bag parameter we take fx = (155.8 ± 1.7) MeV [411 ] . 



B K = 0.725 ± 0.026 

S(x) 

S(x,y) 



4l|. The Inami-Lim functions S(x) and S(x,y) are 42] 
Ax— llx 2 + x 3 3 x 3 \nx 




2(1 



x 



1 



13 1 3 
4 + 21^ - 4(l- 2 /)2 



1 



y-x 



1 3 

4 + 2 1 - a; 



4(1 -x) 



A(l-x)(l-y) 



(35) 



where x = m; l/Myy for all quarks q. 



The predictions for the short-distance QCD factors are: i] c = (1.51 ± 0.24) 43] 



rjd = 0.47 ± 0.04 [44], |45], ^ = 0.58 [46[. The values for r] ct and ?? c have a sizeable 



uncertainty as they are sensitive to the light scale ~ m c where a s is large. The QCD 
correction factor rjt> is given by 



a s (m t ) 



6 / 23 / a s (mb> 



a s (m t ) 



6/21 



Qi s {m t i 
a s {m h > 



6/19 



(36) 



a s (/j,) is the running coupling constant at the scale fi at NLO 48(. Here we assume 
r]tt' = Vt> and 7] ct i = r] ct . 

Using Eqs. ([32]) and ((34]), we get 



-Im 



7^*^)^(20 + 2r ?rf (K^ cs )(^)^(x c ,x i ) 



12 v / 2tt 2 (AM^) cxp 
+ r; t fey ts ) 2 5(x t ) +2^(^V c .)(v5r, l VV.)5(a; e ,^) 

+ 2 Vtt ^%)(V t r d V t/s )S(x t ,x t ,) + r ]t ,(V; d V t , s ) 2 S(x t 



(37) 



The measured value is |e^| = (2.32 ± 0.007) x 10~ 3 [18]. This which mainly puts 
constraints on the combinations V td V ts and V t f d V t ' s , which, to leading order in A, 
depend on A 2 [l — Ce t&ub ] and q 2 [l — (p/q)e^ 5cb '~ 5ub '^}, respectively. 



2. K- 



7T 1 VV 



The flavor-changing neutral-current (FCNC) quark-level transition s — > dvv is 
responsible for the decay K + — > 7i + uu. Unlike other K decays, K + — > -k + vv is dom- 
inated by the short-distance (SD) interactions. The long-distance (LD) contribution 



11 



to Br (K + — > 7i + uu) is about three orders of magnitude smaller than that of the 



SD 



49 



50j. As the decay K + — > -k + vv occurs via loops containing virtual heavy 



particles, it is sensitive to the fourth-generation quark t' . 

The effective Hamiltonian for the decay K + — > -k + vv in the SM4 can be written 

as 

Gp Oi 



eff 



+ V t 1 s V t > d X{x t >) (sd) v _ A (uiUi) v _ A 



(38) 



The function X(x) (x = mf^/M^r), relevant for the t and t' pieces, is given by 

X{x) = r] X X (x) , (39) 



where 



Xn(x) 



X 



2 + x 3x — 6 

+ ~ — mx 



1 — x ' (1 — x) 2 



(40) 



Above, rjx is the NLO QCD correction; its value is estimated to be 0.994 40j. The 
function corresponding to X(x) in the charm sector is X l NL : 



C NL - AB 



(1/2) 
NL i 



(41) 



where C^i and B^^ correspond to the electroweak-penguin and box contributions, 
respectively. The explicit forms of Cnl an d B^^ are given in Refs. 



48 



3- 



The branching ratio of K + — > ti + vv in the SM4 is given by 
Br(K + — > ti + vu) = k 



+ 



A 



A 5 



+ Mv^ vd ) x{xti 



(42) 



where 



3a 2 Br(K + 7r°e+z/) , 

K + = r K+ 7r-n A 



Po(X) 



1 

A 1 



(43) 



We see that Br(K + — > -k + vv) is related to the experimentally well-known quantity 
Br(K + — > TT°e + b>). summarizes the isospin-breaking corrections in relating 
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K + — > -n + vv to K + — > 7r°e + z/; its value is tk+ = 0.901. k + is estimated to be 



(5.36 ±0.026) x HT 11 |52|. 



The measured value is Br(K + — > tt + vv) = (1.7 ±1.1) x 10~ 10 [18|. This which 
mainly puts constraints on the combinations V t * d V ts and V t * d V t ' s , which, to leading 
order in A, depend on A 2 [l — Ce l5ub ] and q 2 [l — (p/q)e l ( Scb '~ s ™ b '* > ], respectively. 

E. The B system 

Here we present various observables in the B system with the addition of a fourth 
generation. 

L B d,s-Bd,s mixin 9 

Meson-antimeson mixing occurs in the SM through the box diagram, and is thus 
sensitive to new heavy particles appearing in the loop. Within the three-generation 
SM, the dominant contribution to B q -B q mixing (q = d, s) comes from the virtual 
top quark. The charm and the mixed top-charm contributions are negligibly small, 
and hence the analysis is simplified considerably. In the SM4, there is an additional 
contribution due to the virtual t' in the box diagram. 

The mass difference AM q is given by AM ~ 2|Mj ? 2 |, where Mf 2 is the virtual 
part of the box diagrams responsible for the mixing. For B®-B® mixing, Mf 2 in the 
SM4 is given by 

f 2 M 2 r ~ ~ ~ ~ 

M? 2 = ^^m B S 6 ji g ^(^^) 2 S(xO±r/H^^ 6 ) 2 ^(x,) 

+ 2 Va iV*V tb )(V; q V t , b )S(x t ,x t ,)] , (44) 

where x t = m 2 /m 2 v , x t > = m 2 /M^ v . The Inami-Lim functions S(x) and S(x,y) 
are given in Eq. (1351) . Here we assume i] t r = i] tt / for simplicity. The numerical 
values of the structure functions S(x t >), S(x t ,x t ') and the QCD correction factor 
r/f for various t' mass are given in Ref. In order to reduce the sizeable non- 
perturbative uncertainties due to the decay constant fs and the bag parameter 
B bq , we consider the ratio AM S / AM d : 



AM, m Bs /2 



A M d m Bd 
13 



C A x M R , (45) 



where 



M, 



R 



Vt(V t * s V tb 


) 2 s(x t ) + r]t/ (v t r s v t , b 


) 2 S(x t ,) 


+ Zrjtt 


(V t * s V tb ) 


(v; s v t , b )s{x t ,x t ,) 


Vt(V t * d V tb 


) 2 s(x t ) + Vt ,(v; d v t , b 


) 2 S(x t ,) 


+ 2r]tt 




{V t 1 d V vh )S{x t ,x t .) 



_ (46) 
with £ = f Ba \fBbsl 'fBjVWd = 1-243 ± 0.028 [4l|. There is less uncertainty in f 



2-3%) than in f B J B bq (~ 7- 



The measured values are 18] 



AM S = (17.77 ± 0.12) ps" 1 , AM d = (0.507 ± 0.005) ps^ 1 , (47) 

whose ratio is sensitive to V t lV ts , V t lV t d, V t f b V t ' s , and V t J b V t fd- These correspond to 
the combinations of the CKM4 parameters A 2 , A 2 (l — Ce lSub ), qre lScb ' , and r(ge ll5ci '' — 
pe i5ub '), respectively, to leading order in A. 

2. CP violation 

CP violation in the quark sector is due to phases in the quark mixing matrix. 



In the three-generation SM, the phase in 
encapsulated in the unitarity triangle 



ormation in the CKM matrix is elegantly 
, whose interior angles are a, /3 and 7. In 



order to test the SM, these angles must be measured in as many ways as possible 
to test for consistency. Unknown strong QCD phases contaminate many of these 
methods; ways of removing these strong phases must be devised in order to cleanly 
measure the weak phases. In the SM4, many of the ways of eliminating the strong 
phases fail, since typically there are multiple amplitudes with different strong and 
weak phases. Here we consider only those constraints from CP observables which 
are free from uncertainties due to the strong phases. 

• Sj/^k s '- The coefficient of sin(AMdt) in the time-dependent indirect CP asym- 
metry in B d — > J/ipKs is given by 

S J/ilKs = 8^20™, (48) 

where 0g* is defined as 

Mf 2 = \M? 2 \e 2 ^ . (49) 
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Thus we have 



good 



In the SM, this is sin 2(3, and is free of strong phases. It is thus a 
observable to constrain the SM4 using Eq. (jHJ). The measured value is 18] 

Sj /i>Ks = 0.672 ± 0.024 , (51) 

which is sensitive to V^ b V t d and V t * b Vt'd, i-e. to the parameter combinations 
A 2 (l — Ce t&ub ) and r{qe lS <= h ' — pe tSub '), respectively. 

Recently, the CDF and DO collaborations measured indirect CP violation in 



B® — > J/ip<f) an d found a 2.2cx deviation from the prediction of the SM 53j. At 
first sight, this seems to indicate a nonzero phase of B a s -B® mixing, and many 
papers have been written exploring the contribution of particular new-physics 



models to this mixing (including the fourth generation J4]). However, there 
could be a significant contribution to this signal from new physics in the decay 



54j. If this is the case, strong phases will play a role. For this reason, the 



constraints from this measurement are not included in the fit. 

• In the SM, 7 = Arg(— V* b V u d) / (V* b V c d) ■ This phase can be probed in tree-level 
decays. By measuring several different decays, it is possible to remove the 
dependence on the strong phase and extract 7. The latest value is [iij 

7 = (75.0 ±22.0)° . (52) 

Because this angle is measured in tree-level decays, its value is unchanged with 
the addition of a fourth generation. Indeed, from Eqs. (jlj) and (JBJ), we see that 

V Vcb v cdJ \ v; b v cd J 

Thus, the phase 5 u b can be constrained through the measurement of the weak 
phase 7, and this observable is included in the fit. 

3. B -»• X s j 

The quark-level transition b — > 57 induces the inclusive radiative decay B — > X s ~f. 
This decay can occur only at the loop level and hence is suppressed within the SM. 
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It has been observed with a branching ratio of (3.55 ± 0.25) x 10 4 55|, in good 

agreement with the NNLO SM prediction of (3.15±0.23) x 10" 4 |56|. Thus, B X s -f 
has a great potential to constrain new-physics models. 

Within the SM, the effective Hamiltonian for the quark-level transition 6—7-57 
can be written as 

AC 8 

n eff = -£vt.vz ci(/0 Qi{n) , (54) 



-^-V ts V t lJ2Ci^)Qi(n), 

" i=l 

where the form of the operators Q 7 ;(/x) and the expressions for calculating the Wilson 
coefficients Cj(/i) are given in Ref. |57|]. The introduction of a fourth generation, in 
addition to the modifications V ts —> V ts and V t b —> V t b, also changes the values of 
the Wilson coefficients C^g via the virtual exchange of the t'-quark. They can be 

written as _ _ 

Vl.V., . 

(55) 



C?l{p) = CrM + ^iCf.(Ai) • 



v; b v ts 

The values of C\ g can be calculated from the expressions for C 7)8 by replacing m t 
by mt'- 

In order to reduce the large uncertainties arising from 6-quark mass, we consider 
the following ratio 

Br{B->X al ) 
Br(B X c eue 



In leading logarithmic approximation this ratio can be written as 40 ] 

2 

Qa\C) ot {m b )\ 2 



R 



VSYu 



v, 



cb 



irf(m c )K(m c 



(56) 



Here the Wilson coefficient C7 is evaluated jit the scale // = m&. The phase space 
factor f(m c ) in BriB — y X c ev) is given by 



f(rh c 



1 — 8m c + &m c — rh c — 24m c In m c 



where m c = m c jm b . K(rn c ) is the 1-loop QCD correction factor 58]: 



K(m c ) 



2a 3 (m fe ) 
3tt 



7T 



T " 



+ 2 



(57) 



(58) 



The values of the branching ratios in R are Br(B — y X s ^) = (3.55 ± 0.25) x 10 4 
Q and Br(B -> X c e^) = 0.1061 ± 0.0016 ± 0.0006 [59]), the ratio being sensitive 
to V t * b V ts and V t J b V t > s , i.e. to A 2 and qre t5cb ' , respectively, to leading order in A. 
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4. B^x s i + r 



The effective Hamiltonian for the quark-level transition b — > s l + I in the SM 
can be written as 

AC 10 

lUft = -7f V ^tl E C &) Qifa) , ( 59 ) 

* i=l 

where the form of the operators Qi and the expressions for calculating the coefficients 



Cj are given in Ref. [57||. The fourth generation, in addition to the modifications 
Vts — > V ts and Vtb —> Vtb, changes the values of the Wilson coefficients (77,8,9,10 via the 
virtual exchange of the if . The Wilson coefficients in the SM4 can then be written 

as _ _ 

CfV) = CM + W^Ct'(fi b ) , (60) 

where % = 7,8,9,10. The new Wilson coefficients C\ (^&) can easily be calculated 
by substituting m t > for m t in the SM expressions involving the t quark. 
The calculation of the differential decay rate gives 

dBr(B -»■ X.l+T) _ a 2 Br(B X c ev) \V* b V ts \ 2 



dz 47r 2 f(m c )K,(m c ) \y cb \ 

where 



2 



[l-zYD(z), (61) 



D(z) = (1 + 2z) (|Cg°T + \Clf\ 2 ) + 4 (l + ~j |C* ot | 2 + 12Re(C* ot C* ot *) . (62) 

Here z = q 2 /m 2 and m g = m q /mb for all quarks g. The expressions for the phase- 
space factor f(rh c ) and the 1-loop QCD correction factor K(m c ) are given Eqs. fl57|) 
and (158]) . respectively. 

The theoretical prediction for the branching ratio of B — > X s l + l~ in the in- 
termediate q 2 region (7 GeV 2 < q 2 < 12 GeV 2 ) is rather uncertain due to the 
nearby charmed resonances. The predictions are relatively cleaner in the low-g 2 
(lGeV 2 < q 2 < 6 GeV 2 ) and the high-g 2 (14.4 GeV 2 < q 2 < ml) regions. Hence we 
consider both low-g 2 and high-g 2 regions in the fit. The branching ratios are 60l. |61| 

Br(B ^ X s l + r ) lowg 2 = (1.60 ± 0.50) x 10~ 6 , (63) 
Br(B -> X s /+ n) highg 2 = (0.44 ± 0.12) x 10~ 6 . (64) 

Both of these branching ratios are sensitive to V t * b Vt s and V t J b Vt> s , i.e. to A 2 and 
qre l5cb ' , respectively, to leading order in A. 
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F. The D system 



In principle, there can be constraints from D°-D° mixing. In the SM, this mixing 
arises due to d, s and b quarks in the box diagram. The b contribution is enhanced 
by a factor of (m| — mjj/(iiij — m^). On the other hand, it suffers a strong CKM 
suppression by a factor of iKbKfelVlKisKisI 2 which is ~ A 8 . Thus, D°-D° mixing 
is dominated by the d- and s-quark contributions. As a result, the mixing is small 
within the SM and hence sensitive to new physics. 

There have been attempts to constrain the CKM4 parameters using D°-D° mixing 



(for example, see Ref. [62] ). However, precisely because the d and s quarks dominate, 
there can be large long-distance (LD) contributions to the mixing. At present, there 
is no definitive estimate of these LD effects. Because of this, we do not have an 
accurate enough prediction for D°-D° mixing, and this measurement cannot be 
incorporated in the fit at present. 

IV. RESULTS OF THE FIT 

We perform the fit to 9 CKM4 parameters, using the observables described in 
the previous section. We define 

Xtotal = XCKM + XlJC + Xzbb + XzAb + xfe K | + Xk+^tt+vD + Xmixing 

+ X 7 + Xb — >X s f Xincl-low Xincl-high J (^5) 

where the exact definition of each \ 2 contribution is given in the Appendix |A] We 
perform the fit at two values of t' mass: m t ' = 400 GeV and m t ' = 600 GeV. In 
addition, we also perform a fit for the 4 parameters of the CKM matrix in the SM, 
in order to check for consistency with the standard fit. The results of these fits are 
summarized in Table [B It may be observed that the \ 2 P er degree of freedom is 
small in each case, indicating that all the fits are good. The goodness of fit does not 
seem to depend much on the masses of the heavy quarks. 



The fit for the SM is consistent with that obtained in Ref. [181 ] ■ As far as the 
parameters of the three-generation CKM matrix are concerned, their best-fit values 
are not affected much by the addition of a fourth generation. However, the allowed 
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Parameter 


SM 


m t i = 400 GeV 


m t , = 600 GeV 


A 


0.227 ±0.001 


0.227 ±0.001 


0.227 ±0.001 


A 


0.808 ± 0.021 


0.801 ± 0.022 


0.801 ± 0.022 


C 


0.38 ± 0.01 


0.42 ± 0.04 


0.42 ± 0.04 




1.16 ±0.06 


1.24 ±0.23 


1.22 ±0.24 


P 




1.45 ± 1.20 


1.35 ± 1.56 


q 




0.16 ±0.12 


0.12 ±0.07 


r 




0.30 ±0.37 


0.19 ±0.27 


^ub' 




1.21 ± 1.59 


1.32 ± 1.76 


<W 




1.10 ± 1.64 


1.25 ± 1.81 


X 2 /d.o.f. 


6.64/14 


6.01/11 


6.06/11 



TABLE I: The results of the fit to the parameters of CKM and CKM4. 

parameter space for C and 5 u b expands by almost a factor of four. This is expected, 
since the constraint on \V U }\ from the 3x3 unitarity is now relaxed. 

On the other hand, the new real parameters p, q, r are consistent with zero, 
which is not surprising since the SM fit is a good one. This also is consistent with 
the observation that no meaningful constraints are obtained on the new phases 5 u y 
and 5 c b'\ since vanishing p, q imply vanishing V u b>,V c b>, respectively, the phases of 
these two CKM4 elements have no significance. 

For m t i = 400 GeV, the maximum values of the parameters (p, q, r) are 
(2.65, 0.28, 0.67) to la. For m t i = 600 GeV, the la upper bounds are 
(2.91, 0.19, 0.46). This indicates that these quantities are indeed 0(1) or smaller, so 
that the expansion in A in the DK parametrization is justified. 

The magnitudes of CKM4 elements are of special interest, since the off-diagonal 
elements are indicative of the mixing between generations. Table [Til gives the al- 
lowed ranges for the magnitudes of CKM4 elements, obtained using the fit results in 
Table HI Clearly the extension to four generations only expands the allowed ranges 
of the CKM parameters, while the allowed values of all of the new parameters of 
CKM4 (except VW) are consistent with zero. 

The combinations of CKM4 matrix elements that control B^-B^ and B s -B s mix- 
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Magnitude 


SM 


ma = 400 GeV 


m t , = 600 GeV 


\V U d\ 


0.9743 ± 0.0002 


0.9743 ± 0.0002 


0.9743 ± 0.0002 


\V US \ 


0.227 ±0.001 


0.227 ±0.001 


0.227 ±0.001 


\V ub \ 


(3.55 ±0.17) x 10~ 3 


(3.90 ± 0.38) x 10~ 3 


(3.91 ±0.39) x 10~ 3 


\V U b'\ 




0.017 ±0.014 


0.016 ±0.018 


\Vcd\ 
1 CU | 


0.227 ±0.001 


0.227 ±0.001 


0.227 ±0.001 


\Vcb\ 


0.9743 ± 0.0002 


0.9743 ± 0.0002 


0.9743 ± 0.0002 


\v cb \ 


0.042 ± 0.001 


0.041 ± 0.001 


0.041 ±0.001 


\Vcb>\ 




(8.4 ±6.2) x 10~ 3 


(6.0 ±3.8) x 10~ 3 


\Vtd\ 


0.0086 ± 0.0003 


0.009 ± 0.002 


0.009 ±0.001 


\Vu\ 


0.041 ±0.001 


0.041 ± 0.001 


0.040 ± 0.001 


\Vtb\ 


1 


0.998 ± 0.006 


0.999 ± 0.003 


\V tb '\ 




0.07 ±0.08 


0.04 ±0.06 


\Vt'd\ 




0.01 ±0.01 


0.01 ±0.02 


\Vt's\ 




0.01 ±0.01 


0.004 ±0.010 


\Vtfb\ 




0.07 ±0.08 


0.04 ± 0.06 


\Vfb'\ 




0.998 ± 0.006 


0.999 ± 0.003 



TABLE II: Magnitudes of the CKM4 elements obtained from the fit. 



ing are V t * b V t d, V t * b V ts , V t >bV t id and VfbVt's- The allowed ranges of these quantities 
are given in Table IIHI It may be observed that here the fourth generation can 
have maximal impact. While |V^Vt s | is little affected, the allowed range of |V^Vfd| 
is increased by up to a factor of 6-7. Moreover, the allowed range of the phase of 
VttVtd is expanded by ~ 10, while that of the phase of V t * b V ts is larger by a factor 
of ~ 20 at mt' = 400 GeV. Since the phases of V t J b V t 'd and V t f b Vt> s are essentially 
unconstrained, they can influence and B s mixing to a large extent. In particular, 
the B s -B s mixing phase can be very large, as suggested by the recent measurements 
from B s — > J/ip(j) decays [63]. The combinations V t * s V t d and V t * s V t >d also contribute 
to rare K decays, and hence significant deviations of these quantities from the SM 
can also leave their imprints in the rare K decays. 
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Quantity 


SM 


m t , = 400 GeV 


m t , = 600 GeV 


\V t tVtd\ 

1 TO LLL l 

^g(V*V td ) 


0.0086 ± 0.0003 
(-21.5 ± 1.0)° 


0.009 ± 0.002 
(-30.4 ± 10.3)° 


0.009 ± 0.001 
(-27.9 ±8.0)° 


\VtlVts\ 

1 TO t " 5 I 

Arg(VSy ts ) 


0.041 ± 0.001 
(-178.86 ±0.06)° 


0.040 ± 0.001 
(-178.12 ± 1.14)° 


0.040 ± 0.001 
(-178.12 ±0.57)° 


Arg(V^V^) 




0.0010 ±0.0015 
(-107.1 ± 106.5)° 


0.0006 ±0.0011 
(-102.5 ± 112.8)° 


Ar g (y; b y t , s ) 




0.0005 ±0.0010 
(37.8 ± 120.3)° 


0.0002 ± 0.0005 
(40.1 ± 174.1)° 



TABLE III: Combinations of CKM4 elements that control mixing in the and B s sectors. 

V. DISCUSSION 

In this paper we consider the extension of the standard model (SM) to four 
generations. Using input from many flavor-physics processes, we perform a \ 2 n t 
to constrain the elements of the 4x4 CKM quark-mixing matrix (CKM4). The fit 
takes into account both experimental errors and theoretical uncertainties. Although 
we do not include the oblique parameters in our fit, we do take values for the masses 
of the fourth-generation quarks that are consistent with the oblique corrections. 

At this stage, several comments are in order. 

• The best-fit values of all three new real parameters of the CKM4 matrix are 
consistent with zero. Since the fit to the SM is also excellent - x 2 /d.o.f. = 
6.64/14, corresponding to a goodness-of-fit of 92% - we must conclude that 
the addition of a fourth generation is not necessary to get a better fit to the 
data. 

• We find V tb = 0.998 ± 0.006 for m t , = 400 GeV and V tb = 0.999 ± 0.003 for 
m t , = 600 GeV. Thus, at 3a, we have V tb > 0.98. Therefore the SM4 cannot 
account for any large deviation of V tb from unity. 

• In many previous analyses, it is mentioned that any mixing between the third 
and fourth generations is small. We find that this is indeed the case - the 
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results of the fit constrain the matrix elements describing the mixing of the 
ordinary and fourth-generation quarks to be |14&'| < 0.06, \Vd>'\ < 0.027, and 
\V tb >\ < 0.31 at 3(T. 

However, the allowed parameter ranges still allow large deviations from the 
SM as far as the magnitudes and phases of the quantities V t * b V tc i and V t * b V ts are 
concerned. With additional new-physics contributions involving VfbVt'd and 
VfbVt's, it may be possible to get significant new-physics signals in B d and/or 
B s mixing, which could be the most incisive probes of the fourth generation. 

The value of \V u b\ required to explain the recent measurement of Br(B + — > 
t + u t ) is 2.8a larger than that obtained from the global fit to \V u b\ otherwise, 



within the SM 



641 ] . Our fit indicates that the best fit for \V u b\ shifts to higher 
values with SM4, and the error on this quantity also increases. As a result, 
the SM4 may be able to account for this measurement much better than the 
SM. 
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Appendix A: The x 2 function 

We define our \ 2 function to be 

Xtotal = XCKM + XUC + Xzbb + XzAb + X\e K \ + Xk+^tt+ vv ' Xmixing 

— >Xs7 Xincl-low Xincl-high " (-^--0 

The components of this function are defined below. 

• For the direct measurements of the magnitudes of the elements, 

2 / |K a |- 0.2255 \2 / |Kd|- 0.97418 \2 / |K fl |-1.04 \2 

XcKM ' ' \ 0.0019 J _ V 0.00027 / * V 0.06 / 

0.011 ) \ 0.00036 J \ 0.0011 J ' ^ ' 



+ 
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/ |Kb'| 2 - 0.0001\ 2 /|y cb /| 2 + 0.136\ 2 
Xuc ~ V 0.0011 / + V 0.125 / 

{\V U \ 2 + \V Vd \ 2 ) + 0.002x2 / (|^| 2 +|^| 2 ) + 0.134 x2 



0.005 I V 0.125 

For the Z — >■ 66 decay, 



0.001 

where i2&b is defined in Eq. ( 12T|) and 



0.216 ^ 

Xzw, = 777^ > ( A4 ) 



0.020 

where A& is defined in Eq. ( )26|) . 
For K mixing, 

2 / Jey|- 0.00232 x2 
*l-' = l 0.00046 J ' (A6) 
where ex is defined in Eq. (|37|) . Here experimental and theoretical errors are 

added in quadrature, 



Next, we have 



2 /Br(K + -> vr+z/z/) - 1.7 x 10- 10 \2 

Xx+^+^p - ^ 1.1 x 10" 10 / ' ^ ' 



with Br(K + ->■ tt+z/z/) as in Eq. ( H2l) . 
In 5-meson mixing, 



2 _ / M fl - 22.20 x2 

Xmixing ^ ) i 

with M# as defined in Eq. (HST) and 

AM S m Bd 1 



For CP violation in Bj — >■ J/%pKs 



— = 22.40 ± 1.04 . (A9) 



2 _ fSj/^Ks ~ 0.672 \ 2 



— , (A10) 
with Sj/^Ks defined as in Eq. ( 150]) . 
For the CKM angle 7 

, ? = ( s ^Mm.y , (A11) 
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• For the radiative decay 



( 



100 R- 0.330 \ 2 
0.041 ) 



(A12) 



with R as defined in Eq. fl56l) . Here experimental and theoretical errors are 
added in quadrature. 

• For the leptonic decay, 



Br(B — » X s Z + /~) loW(? 2 has been obtained by integrating Eq. (I6T1) within the 
limits (1 GeV 2 < q 2 < 6GeV 2 ). Here experimental and theoretical errors are 
added in quadrature. 



where Br{B — > X s l + I )highg 2 has been obtained by integrating Eq. (I6Tj) within 
the limits (14.4 GeV 2 < q 2 < m 2 b ). 
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